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ABSTRACT 

We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar 
population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discov¬ 
ered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate- 
band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these 
LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of £ Wo (Lya) ~ 100-300 A, 

(2) M * ~ 10 ia5 -10 1L1 M e , and (3) relatively low specific star formation rates of SFR/M * ~ 0.03-1 Gyr -1 . 

Three of the 6 MAESTLOs have extended Lya emission with a radius of several kpc although they show very 
compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the 
MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused 
by star formation process including the superwind activity. We suggest that this new class of LAEs, MAEST¬ 
LOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic 
star-formation history. 
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1. INTRODUCTION 

Most of massive galaxies in the present universe are pas¬ 
sively evolvii^galm^^with little on-going star formation 
(e.g., Kauff mann et al.l 120031) . In the current understand¬ 
ing of galaxy evolution, massive galaxies are considered to 
have evolved more rapidly than less massive systems in ear¬ 
lier universe: so-called the downsizing evolution of galax¬ 
ies (ICowie et al.ll 19961) . These massive galaxies have formed 
their stars actively by a cosmic age of a few Gyr (redshift 
z ~ 2 -3), whe n the cosmic star formation rate density peaked 
(e.g. JBouwens et all2011l) . After this epoch, their star forma¬ 
tion stopped and they passively evolved into elliptical galax¬ 
ies seen today. However, the quenching mechanism of star 
formation in these massive galaxies has not yet been under¬ 
stood because the process may have occurred in a relatively 
short time scale, m a king it difficult t o observe such events 
(e.g., Renzini 20091 iPeng et all 1201 Ol iDurkalec et al.l 120151: 
Ma ncini et al.ll2015l) . 

To seek for star-forming galaxies in young universe, the Hy¬ 
drogen hya emission provides the most useful tool. There¬ 
fore, many searches for redshifted hya emission have re¬ 
sulted in the discovery of young galaxies beyond 7, corre¬ 
sponding to an cosm i c age of < 750 Mvr (lOno et al.ll20l2t 
Shibuva et al.l 120121 iFinkelstein et al.l 120 1 3t iKonno et al.l 
2014 Schenker et al] 120141) . Among such Lya emitting 
galaxies (Ly a emitters, hereafter LAEs), those with a very 
large equivalent width (EW), i.e., extremely strong LAEs, 
are particularly important in that they can be galax ies in a 
very early stage o f galaxy formation (e.g., Sch aererl 120031: 
iNagao et aLll2007l) . 

In order to search for them, we have carried out a survey 
for extremely strong LAEs over an unprecedentedly large vol¬ 
ume. While most of the detected objects with strong hya 
seem to be young galaxies with small stellar mass as expected 
for LAEs, we have serendipitously found 6 massive galaxies 
with extremely large EW (Lya) and relatively evolved stellar 
population at z ~ 3, that show no evidence for an active galac¬ 
tic nucleus (AGN). Here we present the physical properties of 
this new population that is expected to be in a transition phase 
between star-forming and passive evolution. In this Letter, 
we use a standard cosmology with Et M = 0.3, Da = 0.7, and 
Hq = 70 km s -1 Mpc -1 . 


2. DATA AND ANALYSIS 

In this study, we use the multi-wavelength da ta set from the 
Cosmic Evolution Survey (COSMOS: IScoville et al. 2007 ). 
Optical imaging data with 12 intermediate-band (hereafter, 
IA-band) filters equipped on Subaru/Suprime-Cam allow us 
to pick up strong emission-line objects by a significant flux 
excess in one of the IA bands. The spectral resolution of 
our IA filters is R = A/A A = 20-26, and the 12 IA filters 
cover the whole optical wavelength range from 4270 A to 
8270 A dTaniguchi et al.l i20l5 ). Therefore, we can search 
for strong LAEs at2.5<z<5.8. Although the details of 
our selection procedure for strong LAEs are given elsewhere 
(M. A. R. Kobayashi et al., in prep.), we briefly summarize it 
in the following. 

At first, from the COSMOS Official Photometric Cata¬ 
log (version 2012, ICapak et aD 120071) . we selected objects 
with a significant (3 a) flux excess in an IA band from the 
frequency-matched continuum estimated by using two ad¬ 
jacent broad-band magnitudes. In order to identify which 
emission line causes the IA-band excess of these objects, 


we applied the public photometric redshift code EAZY 
(iBrammer et al.1120081) to the multi-band photometric data 
from optical to MIR, which incl ude CFHT u* and i* Sub¬ 
aru Bg'Vr'i'z! an d 12 IA bands (|Taniguchi et al.ll2007l) . Ul- 
traVISTA YJHK (IMcCracken et al. 2012) and Spitzer/IRAC 
3.6 //m and 4.5 /xm bands (ISanders et al 1 1120071) . The excess 
IA band and any broad bands whose wavelength coverage is 
overlapped with the excess IA band are excluded from the 
photometric redshift calculation. We adopted a line identi¬ 
fication with the highest probability in the volume-weighted 
redshift likelihood function, and assigned the photometric 
redshift assuming the emission line enters into the effective 
wavelength of the excess IA band. We selected LAEs from 
these strong emission-line objects and then performed the 
spectral energy distribution ( SEP) fitting with the GAL AXE V 
population synthesis model dBruzual & Ch ariot 2003) to esti¬ 
mate the physical properties of the LAEs. In the SED fit¬ 
ting, we assumed the exponentially decaying star formation 
histories with an e-folding timescale of r = 0.01-10 Gyr. 
The Chabrier initial ma ss function fC habrier 2003) and the 
Calzetti extinction law (ICalzetti et alTl2000l) were adopted. 
The excess IA band and any broad bands overlapping with 
the excess band were again excluded. Although other strong 
emission lines such as [O II], [O III], and Ha may enter into 
the JHK bands, we used all JHK -band data in the fitting, be¬ 
cause the effect of such emission lines is not expected to be se¬ 
rious for these bands with wide filter bandpasses. In addition 
to the multi-band photometry used in the photometric redshift 
estimate, we also used the IRAC 5.8 fim and 8.0 /im bands 
to obtain more accurate physical properties such as the stellar 
mass and age. Our survey covers a 1.34 deg 2 area in the COS¬ 
MOS field, that is, the overlapped area between COSMOS 
deep region and UltraVISTA DR1 (IMcCracken et al.ll2012l) . 
The wide survey area and wide wavelength coverage of the 12 
IA bands allow us to search for strong LAEs at2.5<z<5.8 
over a very large volume of 5.5 x 10 7 Mpc 3 . 

As a result, we obtain a sample of 589 LAEs at 2.5 < 
z < 5.3. In this sample, 18 LAEs have both an extremely 
large rest-frame equivalent width of EWo(hya) > 100 A and 
a large stellar mass with M * > 10 10 5 M 0 . Hereafter, we call 
these 18 LAEs MAssive, Extremely STrong hya emitting 
Objects (MAESTLOs). Since our main interest is the star¬ 
forming activity in galaxies, we rejected possible active galac¬ 
tic nucleiJAGN) by using the IRAC color criteria proposed 
bvlDonlev et al] (120121) . We also used the XMM-COSMOS 
dHasinger eta 1 .1120071) . Chandra-COSMOS dHlvis et alJl2009b 
ICivano et al .120 1 2l> and Chandra COS MOS Legacy (F. Civano 
et al., in prep.) and VLA-COSMOS (ISchinnerer et al.ll2007l) 
catalogues to reject AGNs. In total, 12 MAESTLOs turn out 
to show evidence for AGN. Accordingly, we obtain a sample 
of 6 MAESTLOs without evidence for AGb0. Their obser¬ 
vational properties are summarized in Table [TJ 

Their sizes are measured in the excess IA-band (i.e., hya 
imag e) and the COSMOS HST/ACS /p 8 i 4 w-band mosaics 
dKoekemoer et al.ll2007l) . corresponding t o the rest-fram e UV 
continuum, by using the GALFIT code (IPeng et al.ll2002l) . We 
fit the observed surface brightness with an exponential law, 
taking account of the point spread function (PSF) of these 

18 Note that a stacking analysis for the 6 MAESTLOs, corresponding to 
a ~ 650 ksec exposure, results in no detection. The 95% upper limit in the 
0.5-2 keV band is 3.63 x 10 -5 counts s -1 which corresponds to a rest-frame 
luminosity of 7.76 x 10 42 ergs s -1 at z ~ 3. 
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Table 1 

Physical Properties of the 6 MAESTLOs 


No. 

Zphot 

logAf* 

T 

age 

E(B-V) 

EWq (Lya) 

log [SFR/(Mq yr- 1 )] 



(Me) 

(Gyr) 

(Gyr) 

(mag) 

(A) 

Lya 

SED 

1 

3.16 

1111+ 008 

O 99 + 0-08 
^ -0.00 

i-6ii° 0 :f 8 

0 03 +0 ' 06 
u - u:, -o.oo 

240l 2 ° 

1.12 ±0.02 

0 63 +0 ' 26 
u -°- 3 -o.oo 

2 

2.81 

1111+ 004 
AA,11 -0.07 

1 cq+ 1-58 

1 ' J ~— 0.59 

i.w±° 0 f 2 

0 9 Q+ 0-03 
u,z ^-0.02 

306 ±20 

1.20 ±0.02 

i 7Q+0-12 
i>/y - 0.08 

3 

2.81 a 

10 90 +0 00 
iU.yu_ooo 

0 05 +0 0 ° 
u,u:) -o.oo 

n oQ+o.oo 
u - zy -o.oo 

0 19+000 

U,A -0.00 

172 ±6 

1.04 ±0.02 

n qq+0-00 

u -°°-o.oo 

4 

3.24 

10 7 i+ 012 
/A - 0.03 

°- 40 ±o: 4 m 

1 9 Q+ 0-62 

1 - 2 °- 0.14 

0 12 +005 
U,iZ - 0.04 

178 ±16 

0.88 ±0.03 

0 94+0.22 
- 0.16 

5 

2.50 b 

10 54 + 0 08 

iO.J^_ 00 4 

°- 06 iao 2 

0 14+ 0 - 26 
U,A - 0.03 

0 40+ 0 04 
u -^ u - 0-03 

107^“ 

0.68 ±0.03 

1 94 + 0 - 37 

A * y -0.27 

6 

3.16 

10 52 +0 03 

iU.jz _ 0 03 

o-sofSS 

1 90 +00 ° 
1 - yu - 0.10 

0 02 +0 ' 02 
u - uz -0.01 

124 ±16 

0.64 ±0.05 

0 41 + 0.08 
U - 4 A - 0.03 

Note. — The No. is given in order of decreasing the estimated stellar mass. Errors for the quantities correspond to 1 -cr confidence interval (i.e., A % 2 < 1) 
estimated from the SED fitting. In the SED fitting, the templates older than cosmic age at z p hot are n °t used. The entry of 0.00 for these errors indicate that there 


is no parameter grid in A % 2 < 1 around the best-fit model parameter. 
a Zspec = 2.798. 
b Zspec — 2.513. 



2-obs (p-tn) 


Aobs (P m ) 


^obs (prn) 


Figure 1. SEDs and HST and Subaru images of the 6 MAESTLOs. In the bottom main panel, the observed data points are shown by filled circles with error 
bars. The red symbol shows the IA band with a significant flux excess by the Lya emission, while the grey symbols represent the broad bands whose wavelength 
coverage overlaps with the excess IA band. The blue curve shows the best-fit model SED from the GALAXEV library. The data affected by the strong Ly a 
emission (red and grey points) are not used in the SED fitting. The inverted triangles represents the 3cr upper limit for the undetected bands. In the top panels, 
the thumbnails of the HST ACS in the F814W filter, the Subaru Suprime-Cam /'-band image, and the excess IA-band image are shown. Each panel is 12" x 12" 
in size. The ACS images are convolved with a Gaussian kernel with a = 1 pix (= 0"03) to reduce the noise. 
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Table 2 

Size measurements of the 6 MAESTLOs 



ACS F814W 

S-Cam i 


S-Cam IA 

No. 

n il (kpc) 

/unres 

m l (kpc) 

/unres 

m l (kpc) 

/unres 

1 

0.52 ±0.08 

0.060 

< 3.87 b 

0.660 

4.50 ±0.61 

0 

2 

_a 

_a 

4.49 ±2.01 

0.175 

6.68 ±0.70 

0 

3 

1.00 ±0.09 

0 

4.21 ±1.75 

0.115 

7.18±0.80 

0 

4 

0.34 ±0.12 

0.100 

< 3.83 b 

0.750 

< 6.20 b 

0.685 

5 

< 0.57 b 

0.880 

<4.12 b 

0.710 

3.77 ±1.83 

0.050 

6 

< 0.54 b 

0.645 

< 3.87 b 

0.670 

< 5.74 b 

0.410 


Note. — thl is half-light radius and / un res is Errors for rnL are based on 
the Monte Carlo simulation described in the text. 
a Out of the ACS/F814W-band data. 
b Unresolved. 

data. Here we fix the Sersic index to n = 1 because our data 
are not deep enough to resolve the degeneracy between the 
radius and Sersic index for our MAESTLOs. The PSF images 
of the excess IA-band and ACS data are measured by combin¬ 
ing relatively bright isolated stars in each image. Note, how¬ 
ever, that we cannot measure the ACS sizes for MAESTLO 
No. 2 since this object is out of the HST/ACS field. Further¬ 
more, we also measure the sizes of the rest-frame UV con¬ 
tinuum using the /'-band data taken with the same instrument 
as that of the excess IA-band data. Note that the half light 
radius of the PSF in the excess IA bands is 0.75-0.83 arcsec, 
while that in the /' band is 0.51 arcsec. In order to estimate 
the uncertainty in the size measurements including the sys¬ 
tematic effects such as the background fluctuation, we carry 
out the Monte Carlo simulation as in previous studies (e.g., 
IStraatman et alJl2QT5h . After adding the best-fit model pro¬ 
file to the image at 200 random positions around the origi¬ 
nal position (in a 2' x 2' region), we re-measure their sizes. 
The standard deviation of these 200 measurements is adopted 
as the size uncertainty. In order to check whether the object 
is significantly extended or not, we also calculated the frac¬ 
tion of the cases that GALFIT returned the “unresolved” flag, 
/imres, in the 200 measurements. The estimated half-light radii 
and errors of the MAESTLOs together with / unr es are given in 
Table [2] 

We have also carried out additional simulations. In these 
simulations, (1) we convolved model galaxies with the best- 
fit F814W-band light profile with the PSF of the IA-band 
data, (2) then, we added them into the IA-band image, and 
(3) we measured their sizes with GALFIT. We performed 200 
such simulations for the two MAESTLOs with the extended 
Lya emission (Nos. 1 and 3) for which ACS image is avail¬ 
able; note that No. 2 is located out of the ACS coverage. We 
then find that GALFIT returned the “unresolved” flag in most 
cases (185/200 and 169/200 for Nos. 1 and 3, respectively). 
Therefore, we have confirmed that their Lya emission is re¬ 
ally extended in their IA images. In conclusion, the three 
MAESTLOs have an extended Lya emission with a size of 
several kpc. 

3. RESULTS & DISCUSSIONS 

In Figure [T] we show the rest-frame UV-NIR SED of the 
6 MAESTLOs together with their thumbnails in the excess 
I A, /', and ACS /fsi 4 W bands. It is found that they are signif¬ 
icantly bright in the rest-frame NIR wavelengths, leading to 
their large estimated stellar masses of 1 og(M*/M 0 ) = 10.5- 
11.1. Another unexpected property is that they show very red 


rest-frame UV-optical colors despite of their extremely large 
£Wo(Lya); i.e., the MAESTLOs show a relatively strong 
4000 A continuum break in the rest-frame optical as well as 
the Lyman break in the rest-frame far-UV. These continuum 
features allow us to identify the flux excess in a concerned IA 
band as the Lya emission line, resulting in an accurate pho¬ 
tometric redshift for them. In fact, two of the 6 MAESTLOs 
have spectroscopic identifications and their spectroscopic red- 
shifts agree with the photometric redshifts estimated from 
the IA-band excess (Nos. 3 and 5; see Table 1). The strong 
4000 A break observed in the MAESTLOs suggests relatively 
old stellar population in them, and their best-fit stellar ages 
based on SED fitting are 1-2 Gy® Thus these galaxies form 
a completely different population from typical high-redshift 
LAEs with small stellar masses an d young stellar ages (e.g., 
lOno et alJl2010t lHagen et akll2014l) . 

Despite their relatively old stellar population, the MAEST¬ 
LOs have extremely large L Wo (Lya) of ~ 100-300 A. In 
order to compare the star formation rates (SFR) estimated 
from Lya luminosity, SFR{ Lya), with that from SED fit¬ 
ting, SFR(SED), we show SFR(Eya)/SFR(SED) ratios of 
the MAESTLOs as a function of stellar mass in Figure [2] (a). 
Here, we use the lKennicuttl ( 19981) relation between SFR and 
L(Ha) combined with both the L(Lya)/L(Ha) ratio of 8.7 
under the case B recombination and a correction factor con¬ 
verting from the Salpeter IMF into the Chabrier IMF (i.e., 
multiplied by a factor of 0.60). 

For typical LAEs, it is found that the SFR ratio decreases 
with increasing stellar mass (lHagen et al.l2014l) . On the other 
hand, for most MAESTLOs, SFR( Lya) is comparable to 
SFR(SED) and thus their SFR ratios are similar to those of 
typical LAEs with much smaller masses. Therefore, it is sug¬ 
gested that the escape fraction of the Lya emission is rela¬ 
tively high in these galaxies and/or that there are other addi¬ 
tional energy sources besides the photoionization by massive 
OB stars. 

In order to investigate their evolutionary stage, we 
show the distribution of MAESTLOs in the SFR(SED)-M* 
plane together with typical LAEs at z ~ 3 (Nilsso n et all 
20071: iGawiser et all 120071: lLai et all 120081: lOno et alJl201tt 
Hagen et al. 2014) and galaxies at z p h 0 t = 2.5-3.2 in the COS¬ 
MOS field (Figure [2] (b)). Compared to normal star-forming 
galaxies on the main sequence at similar stellar masses and 
redshifts, the MAESTLOs have a smaller specific SFRPI, 
sSFR = SFR/M * ~ 0.03-1 Gyr -1 , suggesting that their star 
formation activities are just ceasing and that they are in a tran¬ 
sition phase from actively star-forming into quiescent galax¬ 
ies. This contrasts with normal LAEs that tend to have a 
sSFR similar with or higher than main-sequence galaxies 
(e.g.. lHagen et al.ll2014l) . 

In order to complementarily investigate the star formation 
histories of the MAESTLOs, we show the rest-frame U — V 
vs. V — J diagram in Figure [3J Comparing the quiescent 
galaxies studied bv lMuzzin et all(l2013l) . we find that our four 
MAESTLOs with a low sSFR (Nos. 1, 3, 4, and 6) are lo- 

19 Note that 2 MAESTLOs (i.e., Nos. 3 and 5) have relatively young ages. 
No. 3 has a very short e-folding timescale and a clear Balmer break in the 
SED, suggesting that its SFR is rapidly decreasing. On the other hand, No. 5 
shows relatively weak Balmer/4000 A break compared to the other MAEST¬ 
LOs and its e-folding timescale is highly uncertain. 

20 One exception is MAESTLO No. 5, which has a weak Balmer break 
as mentioned above. Its sSFR is consistent with the main sequence at the 
redshift. 
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Figure 2. (a) SFR(Lya) / SF R(SED) vs. M* and (b) SFR( SED) vs. M* for 
our MAESTLOs (red circles) and LAEs at z ~ 3 from the literatures (blue 
and gray circles f or the results from stacked LAEs and from individual LAE, 
respectively). For lHagen et all 12014) . we plot only the LAEs at z = 2.5-3.2 
in their sample. Note that the SFRs and M* in the literatures derived with the 
Salpeter IMF are con verted into those with the Chabrier IMF we use. Note 
also that Hagen et a l. (120141) used the SFRs estimated from the dust-corrected 
UV luminosity instead of those based on the SED fitting. In panel (b), we also 
plot the so-called main sequence of star-forming galaxies at z p hot = 2.5-3.2 
in the COSMOS field in gray scale with bin size of 0.1 dex in M* and SFR. 
The scale (the number of galaxies per bin) is shown in the upper-left inset. 
Solid, dot-dashed, and dotted lines show the median, 16 and 84 percentiles, 
and 2.5 and 97.5 percentiles in bins of 0.1 dex in stellar mass. The dashed line 
shows the relation of sSFR = 1 Gyr -1 . The vertical dotted line corresponds 
to log (Af* /Mo) = 10.5, which is one of the criteria of MAESTLO. 

cated around the selection boundary for the quiescent galax¬ 
ies and that their colors are consistent with the model tracks 
where star formation has been recently quenched. There¬ 
fore, this color analysis reinforces our scenario, suggesting 
that they have been recently quenched and are moving into 
the passive evolution phase. Although the colors of the other 
two MAESTLOs are consistent with the star-forming models, 
their colors can also be interpreted as a galaxy that ceased its 
star formation recently. The larger dust contents in these two 
galaxies may be expected if they are in an early phase of the 
superwind activity; i.e., most of dust grains may have not yet 
been blown out by the superwind. We thus infer that MAEST- 



Figure 3. Rest-frame U — V vs. V — J two-color diagram for the MAEST¬ 
LOs. The 6 MAESTLOs are shown by filled circles with their ID numbers. 
The domain of quiescent galaxies studied bv ttMuzzin et al. (204 3) is shown 
with the gray shaded area. Two star formation models (the constant star for¬ 
mation and the single 1 Gyr starburst) are also shown by blue and orange 
dotted curves for Ay = 0 and by blue and orange dashed curves for Ay = 0.5 
mag. 

LOs are in the final stage of massive galaxy formation where 
their SFRs decrease as gas is ejected from the galaxy as the 
superwind. 

As shown in Table [2 the sizes in the rest-frame UV con¬ 
tinuum of MAESTLOs are small (i.e., < 1 kpc). It is note¬ 
worthy that their sizes are very similar to tho se of c ompact 
massi ve quiescent galaxies found at z ~ 2 (Ivan der Wei et al.1 
120141) . implying that the MAESTLOs can be interpreted as 
their progenitors. It has been recently suggested that massive 
compact star-forming galaxies at z ~ 2-3 evolve into com¬ 
pact quiescen t galaxies aft er the cease of their star forma¬ 
tion (e.g. JBarro et alJl20f3 '). Although they are mostly dusty 
galaxies whose sizes are as small as the MAESTLOs, they 
have a younger age of 1.1 Gyr and a higher sSFR of 0.3- 

3 Gyr -1 than MAESTLOs (Barro et al. 2014 ). We therefore 
suggest that they will evolve to passive galaxies through the 
MAESTLO phase. 

Our survey volume (z = 2.42-2.59 and 2.72-3.42) corre¬ 
sponds to 1.4 x 10 7 Mpc 3 and the number density of MAEST¬ 
LOs is 4.3 x 10 -7 Mpc -3 (2.9 x 10 -7 Mpc -3 if we ex¬ 
clude MAESTLOs Nos. 2 and 5 with a relatively high sSFR). 
Thus MAESTLOs may have been missed by previous narrow- 
band surveys because their survey volumes were insuffi¬ 
cient even in a survey with powerful instruments such as 
Subaru/Suprime-Cam fe.g JOuchi et al.ll2008l) . 

Here we compare the number density of MAESTLOs with 
star-forming and quiescent galaxies with M * > 10 10 5 M 0 us¬ 
ing the stellar mass function of galaxies at the same redshift 
range (lllbert et al.112013 ). We then find that the MAESTLOs 
constitute only 0.2-0.6% of star-forming galaxies and ^2% 
of quiescent galaxies. If we assume that all galaxies with 
M * > 10 10 - 5 M 0 pass the phase of MAESTLO when they 
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evolve from star-forming to quiescent galaxies, we obtain the 
duration of this phase, ~ 0.02 x t un { Y (z ~ 3-4) ^ 30-50 Myr 
20-34 Myr if we exclude Nos. 2 and 5), making them a 
rare population. Such a short timescale truncation has been 
recently discusse d based on other observational properties of 
galaxies at z ~ 3 dDurkalec et al.ll2015tlMancini et al.ll2015l) . 

Finally, we mention about the extended nature of Lya emis¬ 
sion. As demonstrated in Section 2, three of the 6 MAEST- 
LOs (Nos. 1, 2, and 3) have extended Lya emission in the 
excess IA band, while all the MAESTLOs are compact in 
the both ACS /p 8 i 4 w-band and Subaru /'-band data (see Ta¬ 
ble [21). The half-light radius of these three MAESTLOs in 
the IA-band data is ^ 4-7 kpc, while that in the ACS data 
is < 1 kpc for all the MAESTLOs (Table [2]). This extended 
nature of Lya emission is intimately related to the observed 
extremely large EWo(Lya). Plausible origins of the extended 
Lya emission is attributed to (1) scattering of Lya photon 
supplied from the central region of each MAESTLO, (2) pho- 
toionized gas by massive star in the central region of each 
MAESTLO, or (3) shock-heated gas driven by a superwind. 
In addition, there are other possible ideas to explain the ob¬ 
served extended Lya emission. The first idea is a projection 
effect by a nearby LAE with a MAESTLO. However, we con¬ 
sider that such a projection effect cannot be the origin of most 
MAESTLOs because of the small number densities of both 
massive galaxies with relatively low sSFR and LAEs with ex¬ 
tremely large EWs at Z r ^ J 3 (e.g., ~ 1.3 x 10 5 Mpc 3 for 
galaxies with M * > 10 105 M 0 and sSFR = 0.03-0.3 Gyr -1 , 
and ^ 5.4 x 10 -6 Mpc -3 for LAEs with EWo (Lya) > 100 A). 
Namely, the projection probability of these objects with simi¬ 
lar redshifts (e.g., Az < 0.5) is expected to be extremely low 
(~ 0.003 such chance alignments within 1" in our survey vol¬ 
ume). Another idea is that a star-forming dwarf galaxy is go¬ 
ing to merge onto these MAESTLOs. In this case, we have to 
explain why a merging dwarf galaxy experiences such active 
star formation. Although we cannot determine which mech¬ 
anism is dominant solely using the present data, future de¬ 
tailed observations of MAESTLOs such as integral field spec¬ 
troscopy will be useful for this issue. 

How the star formation was quenched in high-redshift mas¬ 
sive galaxies is now the most important issue for understand¬ 
ing galaxy formation and evolution. Therefore, large vol¬ 
ume surveys for such massive galaxies with extremely large 
EWoiLya) will become more important in future. 


for their invaluable help, and all members of the COSMOS 
team. We would also like to thank the anonymous referee 
for valuable suggestions and comments. We also thank Alex 
Hagen for kindly providing us with the information of their 
LAEs. This work was financially supported in part by JSPS 
(YT: 15340059, 17253001, 19340046, and 23244031; TN: 
23654068 and 25707010). 


REFERENCES 


Barro, G., Faber, S. M., Perez-Gonzalez, P. G., et al. 2013, ApJ, 765, 104 

Barro, G., Faber, S. M., Perez-Gonzalez, P. G., et al. 2014, ApJ, 791, 52 

Bouwens, R. J., Illingworth, G. D., Labbe, I., et al. 2011, Nature, 469, 504 

Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 

Capak, R, Aussel, H., Ajiki, M., et al. 2007, ApJS, 172, 99 

Chabrier, G. 2003, PASP, 115, 763 

Civano, F., Elvis, M., Brusa, M., et al. 2012, ApJS, 201, 30 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 

Donley, J. F., Koekemoer, A. M., Brusa, M., et al. 2012, ApJ, 748, 142 

Durkalec, A., Fe Fevre, O., de la Torre, S., et al. 2015, A&A, 576, F7 

Elvis, M., Civano, F., Vignali, C., et al. 2009, ApJS, 184, 158 

Finkelstein, S. L., Papovich, C., Dickinson, M., et al. 20l3, Nature, 502, 524 

Gawiser, E., Francke, H., Lai, K., et al. 2007, ApJ, 671, 278 

Hagen, A., Ciardullo, R., Gronwall, C., et al. 2014, ApJ, 786, 59 

Hasinger, G., Cappelluti, N., Brunner, H., et al. 2007, ApJS, 172, 29 

Ilbert, O., McCracken, H. J., Le Fevre, O., et al. 2013, A&A, 556, A55 

Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, MNRAS, 

341, 33 

Kennicutt, R. C., Jr. 1998, ARA&A, 36, 189 

Koekemoer, A. M., Aussel, H., Calzetti, D., et al. 2007, ApJS, 172, 196 
Konno, A., Ouchi, M., Ono, Y., et al. 2014, ApJ, 797, 16 
Lai, K., Huang, J.-S., Fazio, G., et al. 2008, ApJ, 674, 70 
Mancini, C., Renzini, A., Daddi, E., et al. 2015, MNRAS, 450, 763 
McCracken, H. J., Milvang-Jensen, B., Dunlop, J., et al. 2012, A&A, 544, 
A156 

Muzzin, A., Marchesini, D., Stefanon, M., et al. 2013, ApJ, 777, 18 
Nagao, T., Murayama, T., Maiolino, R., et al. 2007, A&A, 468, 877 
Nilsson, K. K., Mpller, P, Moller, O., et al. 2007, A&A, 471, 71 
Ono, Y., Ouchi, M., Mobasher, B., et al. 2012, ApJ, 744, 83 
Ono, Y., Ouchi, M., Shimasaku, K., et al. 2010, MNRAS, 402, 1580 
Ouchi, M., Shimasaku, K., Akiyama, M., et al. 2008, ApJS, 176, 301 
Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 
Peng, Y.-j., Lilly, S. J., Kovac, K., et al. 2010, ApJ, 721, 193 
Renzini, A. 2009, MNRAS, 398, L58 

Sanders, D. B., Salvato, M., Aussel, H., et al. 2007, ApJS, 172, 86 
Schaerer, D. 2003, A&A, 397, 527 

Schenker, M. A., Ellis, R. S., Konidaris, N. R, & Stark, D. P. 2014, ApJ, 

795, 20 

Schinnerer, E., Smolcic, V., Carilli, C. L., et al. 2007, ApJS, 172, 46 
Scoville, N., Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1 
Shibuya, T., Kashikawa, N., Ota, K., et al. 2012, ApJ, 752, 114 
Straatman, C. M. S., Labbe, I., Spider, L. R., et al. 2015, arXiv: 1506.01380 
Taniguchi, Y., Kajisawa, M., Kobayashi, M. A. R., et al. 2015, submitted to 
PASJ 

Taniguchi, Y., Scoville, N., Murayama, T., et al. 2007, ApJS, 172, 9 
van der Wei, A., Franx, M., van Dokkum, P. G., et al. 2014, ApJ, 788, 28 


We would like to thank both the Subaru and HST staff 




